Gravitational collapse of radiation fluid in higher dimensions 
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I. INTRODUCTION 

The problem of gravitational collapse of massive bodies is fascinating and a long standing 
one in general relativity. This problem arises since general relativity predicts vanishing 
pressure gradients against daunting gravitational forces in massive objects (of the order of 
solar masses). Generally, spherically symmetrical object of mass M and radius r related by 
r ~ 2M (in units c = G = 1), undergoes unrestricted collapse [lj. This process may not be 
spherically symmetrical if the collapsing object possesses angular momentum. The theory 
of gravitational collapse has been widely studied. Examples are: Collapse of a massive dust 
cloud [2], charged perfect fluid sphere [3], rotating massive body [HE], r ole of bulk viscosity 
during collapse [6j[7], collapse in the background of cosmological constant [SJ, thin spherical 
shell of dust [9 J , collapse of homogeneous scalar fields [H]] and collapse in higher dimensional 
spacetimes [TT| [12]. 

Almost all the gravitational collapse models lead to the formation of spacetime singular- 
ities generally hidden by one or more horizons. A singularity is a region where invariants 
like Kretschman scalar and curvature scalar diverge. Numerical simulations of gravitational 
collapse of spheroids show that if the collapsing spheroid is sufficiently compact, the sin- 
gularities are hidden inside the event horizon while they become naked (devoid of event 
horizon) if the spheroid is sufficiently large |13J . However, there are some models in which 
the formation of singularity is avoided e.g. if the collapsing star radiates all the matter [14]. 
Another such model is that of a 'regular phantom black hole' which contains Schwarzschild 
like causal structure and the singularity is replaced by the de Sitter infinity [15] . 

There has been huge interest in naked singularities, although their existence is not very 
clear and these are prohibited by the cosmic censorship hypothesis |16j . The existence 
and formation of naked singularities has been suggested for the gravitational collapse in 
self-similar spacetimes [T7j. The visibility of a singularity is possible if there exists a null 
geodesic emanating from the singularity. Then it requires the existence of families of future 
directed non-spacelike curves which emanate from the vicinity of the singularity [IS] . The 
observation of these non-spacelike curves will give sufficient information about the singularity 
itself. Such a mysterious singularity can also be observable if sufficiently strong shearing 
effects near the singularity delay the formation of the event horizon [19]. Another such 
possibility is that of a black hole accreting phantom energy which results in its 'evaporation' 
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and leading to a naked singularity [20H2~i] . Astrophysically, the phenomenon of gravitational 
collapse is manifested in the form of a gamma ray burst in which a super-giant star explodes 
and releases immense heat flux while the stellar core collapses to form a black hole remanent 
[25]. 

The problem of spherical collapse of a null fluid has been studied earlier by Husain [26] 
and recently extended to higher dimensions by Debnath et al |27| . We here investigate 
the same problem using equations of state which are more general than the barotropic EoS 
p = up. These EoS yield interesting behaviors for the evolution of mass of the black hole. 



II. MODELING OF SYSTEM 

We assume an (n + 2)-dimensional spherically symmetric Husain spacetime given by 

ds 2 = - (l - ! ^f^ dv 2 + 2dvdr + r 2 dn 2 n , (1) 

where the radial coordinate is restricted in the range < r < oo and the advanced null 
coordinate v = t — r with — oo < v < oo is called the Eddington coordinate. Here we have 
the energy momentum tensor with two components: null radiation fluid and the matter fluid 
i.e. 

_ rp{n) i rp(m) /r>\ 
± — J-fjIU ' ± {IV 1 \*J 

where 

r#> = <n,i v , (3) 

and 

= (P + P)(^Vu + luVv) + P9v ( 4 ) 

Everywhere in this paper, all Greek indices range from 1 to n + 2. Here Z M = (1, 0, 0, 0), 
and rjfj, = (| (l — -IBr) , —1, 0, 0), are future-like Null vectors, satisfying W A = = 0> 
and l\rj x — — 1. Also a is the is the energy density corresponding to the Vaidya null direction. 
We require the energy momentum tensor to satisfy the energy conditions given by (a) Weak 
and strong energy conditions are: o > 0, p > 0, p > 0. (b) Dominant energy condition 
(DEC) is: a > 0, p > p, p > 0. 
The Einstein field equations are 

Gfx V = T^iv, (5) 
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for the metric (1) with matter field having stress-energy tensor given by 



nm' 



P = 2^> ( 6 ) 



V = -^=T> (7) 

nm , , 

a = . 8 

Here prime ' and overdot . denote differentiation with respect to the parameters r and v 
respectively. For the positive definiteness of p, p and a, we require 

(a) rri > 0, m" < and (b) m > 0. (9) 

We take the following cases of equations of state to solve the field equations (6) - (8): 

1. p(v) = —Ar + Bp(v), 

2. p(v) = (Dr + Er 2 )p(v), 

3. p{v) = Cr k p(v). 

Here A, B, C, D and E are arbitrary constants independent of both r and v. Note that 
the barotropic EoS p = up, represents an asymptotically flat spacetime if the parameter u 
is constrained by ^ < u < 1, while u — 1 represents a charged Vaidya solution [26J. 

A. Case-1: EoS Linear in r 

We consider an EoS which is linear in variable r and is given by 

p(v) = -Ar + Bp(v), (10) 

Using Eqs. (6), (7) and (10), we get 

»/ x . r, nBm'(r.v) , , 

m"(r, v) = 2Ar n (11) 

Solving Eq. (11), we obtain 

„l— Bn 2Ar 2+n 

m(r,v)=C 1 (y) + C 2 (y)- ^ + 77^ ^ (12) 

v ' w K '1-Bn (2 + n)(l + n + Bn) K ' 
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where C\{v) and C 2 {v) are arbitrary functions of time v. Differentiating Eq. (12), w.r.t r 
yields 

9 A r n+1 

m>(r, v) = C 2 {v)r- Bn + r~ _i_ p ; ( 13 ) 

1 + n + im 

while second differentiation gives 

m "(r,u) = -C 2 {v)Bnr- l - Bn + ^ + (14) 

1 + n + Bn 

Also differentiation of Eq. (12) w.r.t v gives 

r l— Bn 

m(r,v) = C 1 (v) + C 2 (v)- — . (15) 

1 — Bn 

Now using Eq. (13) w! > + r- Bn C , 2 (v) > 0. It yields 

C 2 (l+n + Bn)' 



r > 



2A 



l+n+Bn 



(16) 



Now — C2( iln+Bn) — ® ^ either A > and C 2 (l + n + Bn) < or vice- versa. The later 
quantity yields C 2 > and 1 + n + -Bn < and vice versa. Also m" < implies 



2A(n + l)r- „ ___ l _ jU „ 



1 + n + Bn 



Bnr- x - an C 2 {v) < 0, (17) 



which gives 

Further, m > implies 



,i+n+Bn < (l + n + Bn) 

2A(l + n) v ; 



j.l—Bn 



or 



Ci(«) + ^(t;) rr ^>0, (19) 

- > ^ • (20) 

C 2 Bn - 1 v ; 

Now horizon of the metric is obtained by 1 — = 0. It implies m(r\v) = r n ~ l which 

further yields 

„i— Bn 2Ar^ n 

c -w +c 'Wr^ + (2+ „) ( H» +B „) --' < 21 > 

which is an algebraic equation in r. 

In the graphs to follow, we plot the evolution of a black hole mass by considering different 
models (i.e. cases 1 to 3) resulting from the equations of state of a null fluid. Our graphs 
result from different choices of functions Cj, % = 1...6 which are chosen quite arbitrarily 
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including polynomial, trigonometric and the exponential functions, and exhibit various be- 
haviors for the mass parameter m. These functions Cj, then lead to increasing, decreasing 
or fluctuating manner of mass. 

In figure 1, we have chosen C\(v) = e v2 and C^f) = sint> 3 . The constant parameters are 
fixed at A = 3, B = — 8 and n — 4. It is shown that mass increases in steps as r increases 
for v > whereas for v < 0, the mass increases as r decreases. In figure 2, the functions 
are chosen as C\(v) = e" 2 and C 2 (v) = sint> 2 ; while the constants are taken to be A = —3, 
B = — 2 and n = 4. The graph shows that this model is similar to the previous one except 
for the symmetry in v i.e. the mass increases with the increases in radial coordinate r. 



B. Case- 2: EoS quadratic in r 



We now take another EoS which is quadratic in r given by 



p{y) = (Dr + Er 2 )p(v). 



(22) 



The governing equation is 



m"(r, v) 



n(D + Er) m'(r, v). 



(23) 



Solving Eq. (23), we obtain 




(24) 



o 



Also differentiation of Eq. (24) w.r.t r gives 




(25) 



Substitution of (25) in (23) gives 




(26) 



Differentiation of Eq. (24) w.r.t v results 




(27) 



o 
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Now m'{r,v) > =>• C 4 (u) exp(-|nr(2D + £r)) > =>• C 4 (u) > 0. Further, m"{r,v) < 
-n(D + Er)C 4 (v) < 0, when n(D + Er) > 0. Also m(r,v) > implies 



CsW >-expfSVA / e -^. (28) 







C 4 (u) L \ 2E J \ tiE 

Now to calculate the horizon we take m(r,v) = r™" 1 which gives 



(D+Er) 

D 2 n\ I 2 



C 3 (v) + C A (v)e W [^) / e" z cfe = r n ~\ (29) 



o 

„2 



In figure 3, we have chosen C%(v) = e v and C±{v) = cotv . The constant parameters are 
fixed at D = 2, E = 3 and n — 6. Here the mass m possesses symmetry about t> = 0. The 
mass eventually decreases for large v. In figure 4, the functions are chosen as C 3 (v) = e v 
and C±{v) = cscf 2 with the same choice of constants as in Fig. 3. The mass decreases when 
v increases except for the singularity at v — 0. 

C. Case-3: EoS with arbitrary power in r 

Let us now take the EoS 

p(v) = Cr k p(v). (30) 
Here C and k are arbitrary constants. The governing differential equation is given by 

m"(r,v) = -■nJCr k - 1 m'(r,v). (31) 
The solution of the above equation is 

x / \ r„,.(Cnr k \~*^(\ Cnr k \ . . 

-.«) = c.(»)- s Qi(»)(— J r(-,_ ), (32) 



m(r, 



where we have an incomplete Gamma function given by 

1 Cnr k 



k k 



OO 

1 



x*- x e~ x dx. (33) 



Differentiation of Eq. (32) w.r.t r, we have 

m'(r,v) = C 6 (v)exp ■ (34) 



8 



Again differentiating yields 



m"(r,v) = -Cnr k 1 C e (v)exp ^- 



Further differentiation w.r.t v gives 



rh(r,v) = C 5 (v)--C 6 (t) 



Cnr h 



(. 



Cnr k \ 
—J' 



1 Cnr 1 



(35) 



fc' k 



(36) 



Now m'(r,v) > implies Cq(v) > 0. Also m"(r,v) < implies —Cnr k 1 Cq(v) exp( 



CY, 



< 



0, hence Cnr C 6 (i>) > 0. Further m(r, t>) > implies 



C 5 (v) ^ r /' (' '///•' 



r 



C 6 (y) k \ k 
Now horizon of the spacetime is obtained by solving 



1 Cnr k 



C 5 (v) - 



C 6 (v)r 
k 



Cnr k 



fc' fc 



1 Cnr fe 



(37) 



k k 



,n-l 



(38) 



In figure 5, we have chosen C 5 (v) = e" 2 and C 6 (f) = v A . The constant parameters are 
fixed at c = 2, fc = 6 and n = 3. In figure 6, the functions are chosen as C$(v) = e v and 
Cq(v) = f 3 with the same choice of constants as in Fig. 5. The mass increases as v increases 
in both cases. Thus the accretion of the null fluid results in the increase in mass of the black 
hole. 

Also note that the dominant energy condition p > p implies 



r > 



nm 



m" 



which is a general expression. For p = —Ar + Bp, it gives 

2Ar n+1 + cir- i?n (l + n + Bn) 



r > —n 



_2A(n + l)r n - c x Bnr' x ~ Bn {\ + n + Bn) 
For p = (Dr + Er 2 )p(t), the condition Eq. (39) leads to 

1 

r > 



D + Er 

It further gives (D + Er)r > 1 =^> r > 1 and (D + Er)r > 1. The later yields r > , 
In the third case, p = Cr k p{t) the DEC implies 



(39) 



(40) 



(41) 



l-D 
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III. CONCLUSION 

In this paper, we have investigated the gravitational collapse model of higher dimensional 
Husain spacetime. We have obtained three different expressions of mass in the correspond- 
ing three cases. These expressions contain certain functions Cj which needs to be chosen 
arbitrarily since no boundary conditions are imposed on the governing dynamical equations. 
However our choices of these functions lead to some interesting results: In cases 1 and 3, the 
mass of black hole is increasing due to accretion of null fluid. These solutions physically de- 
scribe the inward (ingoing) Husain spacetime [28J. However the solutions obtained in case-2 
describe the outward (outgoing) Husain spacetime since the mass is decreasing. Our solu- 
tions also satisfy the weak and dominant energy conditions which are necessarily satisfied 
in the classical gravity. Moreover, the equations of state chosen here, are generalizations of 
the previously used ones in [26] and hence give a much deeper understanding of the process. 
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FIG. 1: 



The mass m(r, v) is plotted against r and v for the choice of functions C\ and C2 




FIG. 2: 



The mass m(r,v) is plotted against r and v for the choice of functions C\ and C2 



FIG. 3: 



The mass m(r, v) is plotted against r and v for the choice of functions C3 and C4 




FIG. 4: 



The mass m(r, v) is plotted against r and v for the choice of functions C3 and C4 




FIG. 6: The mass m(r, v) is plotted against r and v for the choice of functions C5 and Cq 



